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One of the most challenging open problems in heavy quarkonium physics is the double charm 
production in e~ annihilation at B factories. The measured cross section of e^e~ — > J /xj} + 77c 
is much larger than leading order (LO) theoretical predictions. With the nonrelativistic QCD fac- 
torization formalism, we calculate the next-to-leading order (NLO) QCD correction to this process. 
Taking all one- loop self-energy, triangle, box, and pentagon diagrams into account, and factoring 
the Coulomb-singular term into the cc bound state wave function, we get an ultraviolet and infrared 
finite correction to the cross section of e~ — > J/?/) -I- 77c at -^s — 10.6 GeV. We find that the NLO 
QCD correction can substantially enhance the cross section with a K factor (the ratio of NLO to LO 
) of about 1.8-2.1; hence it greatly reduces the large discrepancy between theory and experiment. 
With rric — 1.4GeV and fi — 2mc, the NLO cross section is estimated to be 18.9 fb, which reaches 
to the lower bound of experiment. 

PACS numbers: 13.66.Bc, 12.38.Bx, 14.40.Gx 
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One of the most challenging open problems in heavy 
quarkonium physics and nonrelativistic QCD (NRQCD) 
is the double charm production in e~^e~ annihilation at 
B factories. The inclusive production cross section of 
J/ip via double cc in e+e^ J/i^cc at ^/s — 10.6GeV 
measured by Belle Collaboration |l| is about a factor of 
5 higher than theoretical predictions including both the 
color-singlet and color-octet Q cc contributions in the 
leading order (LO) NRQCD % Even more seriously, the 
exclusive production cross section of double charmonium 
in e'^ ^ J / Tprjc measured by Belle 0,01 

a[J/^ + Tjc]x B-''-[>2] = (25.6 ±2.8 ±3.4) fb, (1) 

and BaBar0 

a[J/ip + r]c]x B''J''[>2] = (l7.6 ± 2.8l2:i) fb, (2) 

could be larger than theoretical predictions by an order 
of magnitude or at least a factor of 5. Here i?''^[> 2] is 
the branching fraction for the rjc to decay into at least 2 
charged tracks, so Eqs. and Q give the lower bound 
for this cross section. Theoretically, treating charmo- 
nium as a nonrelativistic cc bound state, two indepen- 
dent studies by Braaten and Lee and by Liu, He, 
and Chao 8] showed that at LO in the QCD coupling 
constant and the charm quark relative velocity v the 
cross-section of e+e" J/ipi^c at ^/s = 10.6GeV is about 
3.8 ^ 5.5fb (depending on the used parameters, e.g., the 
long-distance matrix element, rric and a^). In compar- 
ison with Eq. or Eq. JJJ, such a large discrepancy 
between theory and data may present a challenge to our 
current understanding of charmonium production based 
on NRQCD and perturbative QCD. 

Some theoretical studies have been suggested in order 
to resolve this large discrepancy problem. In particular, 
Bodwin, Braaten, and Lee proposed 0,01 that processes 



proceeding via two virtual photons may be important, 
and Belle data for J/ip + rjc might essentially include 
the J/tjj + J 1^1) events which were produced via two pho- 
tons. Brodsky, Goldhaber, and Lee suggested that since 
the dominant mechanism for charmonium production in 
e"'"e~ annihilation is expected to be the color-singlet pro- 
cess e"''e~ — > ccgg, the final states observed by Belle 
might contain J/t/; and a M ~ 3GeV spin- J glueball 
Qj (J = 0,2) [Ij. Motivated b y th ese proposals. Belle 
presented an updated analysis and ruled out the 

J/ip + J/'4' and spin-0 glueball scenarios. Ma and Si stud- 
ied this process by treating the charm quark as a light 
quark and using light-cone distribution amplitudes to pa- 
rameterize nonperturbative effects related to the inner 
structure of charmonium |l3|. Similar approaches were 
also considered by Bondar and Chernyak |1J|. But the 
enhanced cross section is sensitive to the specific form of 
quark distributions. Hagiwara, Kou and Qiao obtained a 
result consistent with ReLjTj and Ref. 8j, and conjectured 
that higher-order corrections in as may be huge |l5l| . 
There are also other suggestions to resolve the double 
charmonium problem, and a comprehensive review on 
related topics and recent developments in quarkonium 
physics can be found in Ref. 

In order to further clarify this problem, in this pa- 
per we present a result for the next to leading or- 
der (NLO) QCD correction to the process of e+ + 
e~ ^ J/ip + rjc- As is known, the NLO QCD correc- 
tions are important for quarkonium production in in- 
elastic J/ip photoproduction 17], in Jf-ip plus jet and 
plus prompt photon associated production in two pho- 
ton collisions (1 81] ■ and in gluon fragmentation functions 
for heavy auarkonium|l9l| . 

At LO in tts, J/i^ + Vc can be produced at order a^a^, 
for which we refer to e.g. RefQ. There are four Feynman 
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j/i, 




+2 flipped diagrams 



FIG. 1: Bom diagrams for e (fci)e^(A:2) J / ^p{2pi)rjc{2p2) ■ 



diagrams, two of which are sliown in Fig.^ and the other 
two can be obtained by reversing the arrows on the quark 
hnes. Momenta for the involved particles are assigned as 
e-(fci)e+(fc2) ^ J/V^Zpi) +77c(2p2)- Using the NRQCD 
factorization formalism, we can write down the scatter- 
ing amplitude in the nonrelativistic limit to describe the 
creation of two color-singlet cc pairs at short distances, 
which subsequently hadronize into J/tp + rjc at long dis- 
tances in the e"*"e~ annihilation process. (Note that here 
the color-octet cc contribution is of higher order in v and 
therefore negligible). Choosing the Feynman gauge, we 
get the amplitude of Born diagrams 
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^2)^, Be — I is the electric charge of 



where s = (ki 
the charm quark, p is the Lorentz indices of the virtual 
photon, £ is the polarization vector of J/ip- 2pi and 2p2 
are the momenta of J/ijj and rjc respectively. i?s(0) is 
the radial wave function at the origin of the ground state 
charmonium J/^ and rjc- 

At NLO in a^, the cross section is 

do- OC \MBorn+ MnLO^ 
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The self-energy and triangle diagrams all correspond to 
propagators and vertexes of Born diagrams. There re- 
main twenty-four box and pentagon diagrams. Twelve 
diagrams of them are shown in Fig. |21 The upper cc 
hadronize to J/tp, and the lower to rjc- The other twelve 
diagrams are obtained by reversing the arrows on the 
quark lines. Specially, the associated diagram with Pen- 
tagon N12 exists only by reversing the arrows on the 
lower quark lines which hadronize to rjc- 

The self-energy and triangle diagrams are in general 
ultraviolet (UV) divergent; while the triangle, box, and 
pentagon diagrams are in general infrared (IR) divergent. 
Box N5 and N8 and Pentagon NIO, which have a virtual 
gluon line connected with the cc in a meson, also con- 
tain the Coulomb singularities due to the exchange of 
longitudinal gluons between c and c. In the practical 
calculation, the IR and UV singularities are regularized 
with D = A—2e space-time dimension, and the Coulomb 
singularities are regularized by a small relative velocity v 



between c and c [iTjj v = |Pic — Pi^l/™ , defined in the 
meson rest frame. For the Coulomb-singular part of the 
virtual cross section, we find 



|i?s(0)|V" 




In the second step, the Coulomb-singularity term has 
to be factored out and mapped into the wave functions 
of J/ip and 77c. For the LO expressions of operators 



Qj/lp 



and (O'^- 



.(1) 



are associated with 



i?s(0), and the NLO are proportional to nasCp/v j^. 
And the two operators give a factor of 2 at 0{as) , re- 
sulting in just the Coulomb-singular term in Eq. (jSJ. 

The self-energy and triangle diagrams contain UV sin- 
gularities, which are removed by the renormalization of 
the QCD coupling constant , the charm-quark mass m 
and field tp, and the gluon field A^. Similar to the renor- 
malization scheme in Ref.|lJ|(see also we define 

Zggs, Z„m, tP°= ^^j, A^^ ^/Z^A^, (6) 
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where the superscript labels bare quantities and Zi = 
1+SZi, with i = g,m, 2, 3, are renormalization constants. 
The quantities dZi are oiO{as) and they contain UV sin- 
gularities and finite pieces which depend on the choice of 
renormalization scheme. We define Z2 and Z^ in the on- 
mass-shell (OS) scheme, and Z3 and Zg in the modified 
minimal-subtraction (MS) scheme 

2 , 47r/^2 
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where fi is the renormalization scale, je is the Euler's 
constant and /3o = (11/3)Ca — (4/3)TFn/ is the one- 
loop coefiicient of the QCD beta function, and Uf is the 
number of active quark fiavors. There are three massless 
light quarks u,d,s and one heavy quark c, so n/ — A. 
Color factors are given by Tp = 1/2, Cp = 4/3, Ca = 3 in 
SU{3)c- Differing from Ref.[li|, we take the MS scheme 
for Z3 with no external gluon legs and set nj — A. In 
this scheme, we do not need to calculate the self-energy 
on external quark legs. It turned out that the difference 
for the calculated cross section in different schemes is of 
order of next to next to leading order and can therefore 
be neglected in the NLO result. In the NLO corrections 
we should use the two-loop formula for as{p.), 

Qs(m) _ _J_ phhiL^ 

An ~ " (]IL^ ' 



(8) 
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FIG. 2: Twelve of the twenty-four box and pentagon diagrams for e~(fci)e^(A;2) J / ^p{2pi)rjc{2p2) ■ 



where L = ln(^VA^cD)> and /3i = (34/3)Ca^ - 
ACpTrnf — {20/3)CATFnf is the two-loop coefficient of 
the QCD beta function. 

Pentagon diagrams Nil and N12 can be reduced to in- 
tegrals with a lower number of external legs directly, since 
there are only two independent momenta. Then they can 



be calculated the same way as box diagrams. To treat 
Pentagon NIO in Fig. |21 we need to calculate the five- 
point function EQ[pi,2pi,—p2,—2p2,m,0,m,0,m\, and 
the finite term eI™ ^ where 



Eq — E, 



-Do[-pi,-pi -p2,Pi,0, m,0,TO] + -Do[p2,pi +P2,-P2,0,TO,0, m], 



(9) 
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^o[pi+P2, Pi+2p2, -Pi, 0, 0, TO, to] 



d^g 



2/s(s/2 - Aq ■ pi + Aq ■ p2 - Sm^ 



{2Tr)D (g2„m2)(g-fpi)2((q+2pi)2-TO2)(9-p2)2((9-2p2)2-TO2) 
2V4;^?^tan-i^^- Viln(;;^) ^ 2iAm^-s)yhan--l^+^s {zni3m^-s)+{s-4rn^)H^)) 

Sim^TT^ {Am^ — s)s5/2(16to2 — s)~^ ' 



(10) 



where the IR- and Coulomb-finite term iJg™ is calcu- 
lated with dimension Z? = 4 and velocity v — 0, and 
ln(— s/m^) = ln(-(s -|- iO)/TO^) = ln(s/TO^) — in 
"Pi 



In 



Eq. © the -DohPi, 
4 



Do -Co[-pi,pi,0,TO, to] + — 



P2,Pi, 0, m, 0, to] term is, 
2iTT ~ 21n4 



(11) 



This term will appear in BoxN5, N8. The other IR- 
divergence terms can be calculated like that. Then 
all the IR-divergence terms become Cq[pi, ~p2,0,m,m] 
and Co[pic, — Pic, 0, TO, to]. We find that BoxN3,N6, N7 
and PentagonN12 are IR-finite respectively, and sum of 
BoxNl + N2 + N4 + N9 is IR-finitc, and IR-divergence 
term of Pentagon Nil is canceled by vertex dia- 
grams. IR-divergence and Coulomb-singular terms of 
BoxN5-|-N8 and Pentagon NIO are all related to the 
Co[pic, -pic, 0,to,to] term. With v \pu:-pis\/m — > 0, 



Co = 



f Airy? 
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IR-divergence terms of Box N5 
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Pentagon NIO are 



canceled by counter terms, and the Coulomb singularity 
is mapped into Rs{0). UV term is canceled by counter 
terms. Then the final NLO result for the cross section is 
UV-, IR-, and Coulomb-finite. Details of the calculation 
can be found in a forthcoming paper. 

We now turn into numerical calculations for the cross 
section of e+ + — > J/tp + Vc- To be consistent 
with the NLO result the value of the wave function 
squared at the origin should be extracted from the lep- 
tonic width at NLO of (see e.g. '4]): |i?s(0)p = 
[(9TO2/^)/(16a2(l - ACFas/TT))]T{J/i;^^ e+e"). Using 
the experimental value 5.40 ± 0.15 ± 0.07 KeV 113, we 
obtain |i?s(0)P = 0.978GeV^, which is a factor of 1.21 
larger than 0.810GeV'^ that was used in Refs.|3,ll] from 
potential model calculations. Taking toj/^ = to^^ = 2to 

(in the nonrclativistic limit), to — 1.5 GeV, A— — 

^ 'MS 

338MeV, with Eq.® we find aJu) = 0.259 for ^ = 2m 
(these are the same as in Ref.|a|, except here a larger 
|i?s(0)P is used), and get the cross section in NLO 



(7(e^ 



J/ij + Vc) = 15.7fb, 



(13) 
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m = 1.4GeV and /i = 2to, the NLO cross section is es- 




timated to be 18.9 fb, which reaches to the lower bound 
of experiment. 
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FIG. 3: Cross sections as functions of the renormalization 
scale fi. Here \RsiO)\'^ = 0.978GeV^ A = 0.338GeV, ^ = 
lO.GGeV; NLO results are represented by solid lines and LO 
one by dashed lines; the upper line is for m = 1.4GeV and 
the corresponding lower line is for m = 1.5GeV; the upper 
straight line denotes the central value measured by Belle in 
Eq.Q and the lower straight line by BaBar in Eq.l|5J. 

which is a factor of 1.96 larger than the LO cross section 
8.0 fb. If we set fj. = m and = \/s/2, then as = 
0.369 and 0.211, which result in the cross section 27.5 
and 11.2 fb respectively. If we set m = 1.4GeV and 
fi — 2m, the cross section is 18.9 fb and 9.2 fb at NLO 
and LO respectively. (Our LO result is also consistent 
with Rcf. t] if we take their smaller value for |i?5(0) p and 
/i = y/s/2.) In Fig. 13 we show the calculated e+ + ^ 
J/tp + rjc cross sections at LO and NLO as functions of 
the renormalization scale /i with two mass values rric — 
1.4 GeV and 1.5 GeV, as compared with the Belle and 
BaBar data. We see the NLO QCD correction enhances 
the cross section by about a factor of 2, despite of existing 
theoretical uncertainties. 

The relativistic corrections may further significantly 
enhance the cross section [ilj (see also 0)- The reason 
for the enhancement is quite obvious that in Fig. ^ the 
virtuality of the giuon takes its maximum value of = 
s/4 in the nonrelativistic limit, and taking account of 
the relative momentum between the charm quarks in the 
charmonium will lower the value of the gluon virtuality. 

In conclusion, we find that by taking all one-loop self 
energy, triangle, box, and pentagon diagrams into ac- 
count, and factoring the Coulomb singular term associ- 
ated with the exchange of longitudinal gluons between c 
and c into the cc bound state wave function, we get an 
ultraviolet (UV) and infrared (IR) finite correction to the 
cross section of e+e^ J/ip + rjc at y/s = 10.6 GeV, and 
that the NLO QCD correction can substantially enhance 
the cross section with a K factor (the ratio of NLO to 
LO ) of about 1.8-2.1; and hence it crucially reduces the 
large discrepancy between theory and experiment. With 



[1] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 89, 
142001 (2002). 

[2] P. Cho and A.K. Leibovich, Phys. Rev. D 54, 6690 
(1996); F.Yuan, C.F. Qiao, and K.T. Chao, Phys. Rev. 
D 56, 321 (1997); tbtd, 1663 (1997); S. Baek, P. Ko, 
J. Lee, and H.S. Song, J. Korean Phys. Soc. 33, 97 
(1998); V.V. Kiselev et ai, Phys. Lett. B332, 411 (1994); 
K.Y.Liu, Z.G.He, and K.T.Chao, Phys. Rev. D 68, 
031501(R)(2003). 

[3] K.Y. Liu, Z.G. He and K.T. Chao, Phys. Rev. D69, 
094027 (2004). 

[4] G. T. Bodwin, E. Braaten, and G. P. Lepage, Phys. Rev. 

D 51, 1125 (1995); 55, 5853(E) (1997). 
[5] P. Pakhlov, hep-ex/0412041 

[6] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 72, 

031101 (2005) . 
[7] E. Braaten and J. Lee, Phys. Rev. D 67, 054007 (2003); 

D 72, 099901(E) (2005). 
[8] K. Y. Liu, Z. G. He and K. T. Chao, Phys. Lett. B 557, 

45 (2003). 

[9] G. T. Bodwin, J. Lee and E. Braaten, Phys. Rev. Lett. 
90, 162001 (2003). 
[10] G. T. Bodwin, J. Lee, and E. Braaten, Phys. Rev. D 67, 
054023 (2003). 

[11] S. J. Brodsky, A. S. Goldhaber and J. Lee, Phys. Rev. 

Lett. 91, 112001 (2003). 
[12] K. Abe et a/., (Belle Collaboration), Phys.Rev. D70 

(2004) 071102. 

[13] J.P Ma and Z.G. Si, Phys.Rev. D70 (2004) 074007. 
[14] A.E. Bondar and V.L. Chernyak, Phys. Lett. B612, 215 

(2005) . 

[15] K. Hagiwara, E. Kou and C. F. Qiao, Phys. Lett. B 570, 
39 (2003). 

[16] N. Brambilla et al, hep-ph/0412158 

[17] M. Kramer, Nucl. Phys. B 459 (1996) 3. 

[18] M. Klasen, B.A. Kniehl, L.N. Mihaila, and M. Stein- 

hauser, Nucl. Phys. B713 (2005) 487; Phys. Rev. D71 

(2005) 014016. 

[19] E. Braaten and J. Lee, Nucl. Phys. B586 (2000) 427. 
[20] S.Eidelman, et al. [Particle Data Group], Phys. Lett. B 
592 (2004) 1. 

[21] Z.G. He, Y. Fan, and K.T. Chao, to be submitted. 



